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ABSTRACT 

New Chandra X-ray and Herschel FIR observations enable a multiwavelength study of ac- 
tive galactic nucleus (AGN) heating and intracluster medium (ICM) cooling in the brightest 
cluster galaxy (BCG) of Abell 2597 (z = 0.0821). The new Chandra observations reveal the 
central <; 30 kiloparsec X-ray cavity network to be more extensive than previously thought, 
and associated with enough enthalpy to theoretically inhibit the inferred classical cooling 
flow. Nevertheless, we present new evidence, consistent with previous results, that a moder- 
ately strong residual cooling flow is persisting at 4%-8% of the classically predicted rates in a 
spatially structured manner amid the feedback-driven excavation of the X-ray cavity network. 
New Herschel observations are used to estimate warm and cold dust masses, a lower-limit 
gas-to-dust ratio, and a star formation rate consistent with previous measurements. [O l] and 
CO(2-l) line profiles are used to constrain the kinematics of the ~ 10^ Mq reservoir of cold 
molecular gas. The cooling time profile of the ambient X-ray atmosphere is used to map 
the locations of the observational star formation entropy threshold as well as the theoretical 
thermal instability threshold. Both lie just outside the < 30 kpc central region permeated by 
X-ray cavities, and star formation as well as ionized and molecular gas lie interior to both. 
The young stars are distributed in an elongated region that is aligned with the radio lobes, and 
their estimated ages are both younger and older than the X-ray cavity network, suggesting 
both jet- triggered as well as persistent star formation over the current AGN feedback episode. 
Bright X-ray knots that are coincident with extended Lya and FUV continuum filaments mo- 
tivate a discussion of structured cooling from the ambient hot atmosphere along a projected 
axis that is perpendicular to X-ray cavity and radio axis. We conclude that the cooling ICM is 
the dominant contributor of the cold gas reservoir fueling star formation and AGN activity in 
file Abell 2597 BCG. 

Key words: galaxies: clusters: individual: Abell 2597 - galaxies: active - galaxies: star for- 
mation - galaxies: clusters: intracluster medium - galaxies: clusters: general 
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1 INTRODUCTION 

For a subset of galaxy clusters with sharply peaked X-ray surface 
bright ness profiles, the intracluster medium (ICM, e.g., |Sarazin| 
1986 1 can cool via bremsstrahlung processes from > 10' K to 
<^ 10** K on timescales much shorter than a Gyr within a radius 
of ~100 kpc. Simple models predict that runaway entropy loss by 
gas within this radius accompanies subsonic, nearly isobaric com- 
pression by the ambient hot reservoir, driving a long-lived classi- 
cal cooling flow onto the central brightest cluster galaxy (hereafter 
BCG, see cooling flow reviews by Fabian 1994, Peterson & Fa bian| 
|2006j ). In these "cool core" (CC) clusters, catastrophic condensa- 
tion of the ICM should drive extreme star formation rates (10" - 10^^ 
Mq yr"') amid massive cold gas reservoirs (~ lO'^ Mq) in the 
BCG, and high resolution X-ray spectroscopy should detect bright 
coolant lines stemming from <; lO' Mq cascades of multiphase gas 
condensing from the hot atmosphere. Yet three decades of searches 
for these expected cooling flow mass sinks returned with results 
that were orders of magnitude below predictions, consistent only 
with residual cooling at ~ 1 - 10% of the expected rates (e.g., [Allen 
[De Young & Rob erts||1974{ |Peterson ]pi978; Haynes eTaT 
fen et al. 1978, Shostak et al.'1980','O'Dea & Baum igsvf 
McNamara & O'Connell 1989, O'Dea et al. 1994b, AUen 1995] 
O'Dea & Baum 1997, O'Dea et al. 1998; Peres et al. 1998; Mittaz] 
etal. 2001, Tamura et al. 2001 , Peterson et al. 2001 , Xu et al. 2002) 



T959 



1978 



Sakelhou et al.|2002 ; Edge & Frayer 2003; Peters on et al.|2003| 



Bregman & Lloyd-DavieS|2006(|Sanders et al.|2008^ . 



Attempts to reconcile these results with the high X-ray lumi- 
nosities and long lifetimes associated with the CC phase have in- 
voked non-gravitational heating mechanisms to inhibit or replenish 
an average ~ 90% of ICM radiative losses over the cluster lifetime 
(e.g., review by [Peterson & Fabianf 2006i. One promising candi- 
date for quenching cooling flows in the innermost regions of cool 
cores is a feedback loop regulated by the mechanical dissipation 



of active galactic nucleus (AGN) power (e.g., Rosner & Tucker 
1989; Baum & O'Dea 19"9T] [Churazov et al.||2002| IBirzan et~aL 



2004; Rafferty et al. 2006; Dunn & Fabian''2006','Best et al.|2007| 
Edwards et al..2007,,Mittal et al. 2009 reviews by McNamara&j 
2007 2012 pabianl2012J. The paradigm is supported by 



strong circumstantial evidence, including observations of kilopar- 
sec scale X-ray cavities in spatial correlation with radio emission 
associated with AGN outflows (Bohringer et al. 1993; Fabian e t al.| 
20001 '2006', 'Churazov e t aL]|2001, ,McNamara et al.,.2000, 200r 



Blanton et al. 2001 ; Nuls en et al.[2005l|Ftarm"an et al.|2005l 



2007 



Birzan et al.n2004^ . The power of these radio sources are statis- 



tically correlated with X-ray luminosity from within the cooling 
radius, and the most radio-loud BCGs are ubiquitously associated 
with the strongest cool cores (e.g., |Burns[1990[|Birzan et al.|2004[ 



Rafferty et aL|[2006l pitFal et al.||2009| [Sun|[2009[ [Hudson et al.] 



2oTol r 



Observationally supported simulations show that the outflow- 
ing plasma can drive sound waves and subsonically excavate cavi- 
ties in the thermal gas, which then buoyantly rise, entrain magnetic 
fields and colder ISM phases, and locally thermalize enthalpy asso- 
ciated with their inflation as ambient gas moves to refill their wakes 



(e.g., 'Begelman"2001 Ruszkowski & Begel man|[2002[ Churazov| 



^al. 2002; Reynold set al. |2002[ [Fa bian 2003"| Fabian'et^ 
IsTrzan et al. 2004 , Robin son et al.|20 04 , Dursi & Pfrommer 



2008 



2008 



[Gitti et aT.,201 1 ). The phenomenon is likely episodic at a rate cou- 
pled to the AGN duty cycle, and total energy input when summed 
over the cluster lifetime can range from ~ 10'^ -10" ergs. In prin- 
ciple, this is enough to replenish the total energy budget of ICM ra- 



diative losses for an average of the CC cluster population, although 
the spatial distribution and thermalisation of this energy is one of 
several important problems that challenge the model (e.g., |McNa-| 
[mara & Nulsen|2007| [2012 ^, and references therein). The physics 
coupling AGN mechanical energy to ICM entropy remain poorly 
understood, and thermal conduction, gas sloshing, and dynamical 
friction likely play additional important roles (e.g., 'Sparks et al. 
1989, Ruszkowski & Begelman 2002, Voigt et al. 2002, Brighenti 
& Mathews 2003, El-Zant et al. 2004; Soker et al. 2004, Voigt & 
Fabian 2004, Parrish et al. 2010, Morsony et al. 2010„ZuHone et al. 



2010| Blanton et al.[201 1[ Sparks et al.[2009 2012 1. 



1.1 Residual cooling and star formation amid AGN feedback 

Critical details of the heating and cooling feedback loop are en- 
coded in the mass, energy, and timescale budgets of the low tem- 
perature, high density gas phases preferentially found in CC BCGs, 
such as filamentary forbidden and Balmer emission line nebulae 
and ~ 10' -10'° Mq rese rvoirs of both v ibrationally e xcited and 
molecular gas (e.g., Hu et ak |l985 



colcf] 

etaLT 



Baum 



198^ 



Heckman 



i989';'Jaffe & Bremer' 1997', 'Donahue et al. 2000; Edge'2001 
Edge et al. 2002, Jaffe et al. 2001 ; Wilman et al. 2002 , McNamara] 
et al. 2004, Jaffe et al. 2005 , Rafferty et al. 2006, Egami et al. 2006[ 
Edge et al..2010b ,a ^ Oonk et al.|2010,,Mittal et al.,2011,, Wilman] 



et al.|2009|[20TT] [Salome et aL|2006|]20TTl|Lim et al.|2012^ . 

While residual ICM condensation can contribute a major frac- 



tion of the mass budget for these phenomena (e.g., Baum|]1987 
Heckman et al. 1989, Cavagnolo et al. 2008; Rafferty et al. 2008 
iQuillen et al.|2008i|O'Deaet al.,2008„Hudson et al.,2010„McDon 



[aid et al.|2010 2011b[a) , their temperatures and ionization states 
are often inconsistent with the 10 ^ T 10"* K phases of a purely 
radiative cooling flow (e.g., [Donahue & Voit|]l991[ [Voit & Don-[ 
[ahue[[1997] l. Thermal conduction and suprathermal electron heat- 
ing now appears to be important (IFerland et al. '2008' 12009 [Don-[ 



ahue et al.|20"TTl[Fabi an et al. 201 1 , Oonk 201 1 , Mittal et aL 



20111 



Sparks et al. 2009, 2012, Johnstone et al. 2012, in press), and a 
critical role is played by the clumpy and filamentary distributions 
of star formation ongoing amid these cold reservoirs on <^ 30 kpc 
scales (e.g., 'Johnstone et al.'1987';'Romanishin'1987"McNamara & 
O'Connell 1989; Hu 1992; Crawford & Fabian 1993; McNamara 
et al. 2004 , Hansen et al. 1995 , Allen 1995 , Voit & Donahue 1997 
Smith et al ._1997,,Cardiel et al..l998,,Hutchings & Balogh.2005 



Oegerle et al.|2001[ [Mittaz et ai.|200l[ 



O'Dea et 



al.[[2001 



2004 



McNamara et al.''2004'; 'Hicks & Mushotzky'2005'; 'Rafferty et al 
2006; Bildfell et al. 2008 , Raffer ty et al. 2008, O'Dea et al. 2008[ 
Voit et al. 2008 ; Loubser et al. 2009| [O'Dea et al.[2010|[Oonk et al' 

5oii 



.. .McDonald et al..2010„ 201 lb 



Estimated star formation rates range from a few to tens of so- 
lar masses per year, and strongly correlate with upper limits on 
spectroscopically derived ICM mass deposition rates, as well as 
CO-inferred molecular gas masses (e.g., [O'Dea et al. 2008). This 
suggests a direct causal connection between reduced cooling flows 
and star formation, as BCGs with young stellar populations are al- 
ways found in cool core clusters ( [Bildfell et al.[2008|[Loubser et al.[ 



' Throughout this paper, we will use "hot" to describe 10^ < T < 10** K 
(X-ray bright) ICM/ISM phases, "warm" to describe 10* <; T < 10^ K (op- 
tical and UV bright) components, and "cold" to describe 10<^r<^10''K 
(N/M/FIR bright) components. We will not discuss the critically important 
10^ — 10^ K regime at great length in this paper. 
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Table 1. A summary of the new and archival observations used (directly or referentially) in this analysis. This table also applies to the analysis presented in 
[Tremblay et al.|2012a| (1) Facility name; (2) instrument used for observation; (3) configuration of instrument / facility; (4) wavelength regime / spectral line 
observed or filter used; (5) exposure time; (6) facility-specific observation or program ID; (7) date of observation; (8) principal investigator (PI) of observation 
or program; (9) reference to publication where the data were first published — 1 : [McNamara et al. (2001 1; 2: Oonk et al. (2011i; 3:'0'Dea et ar](2004) ; 4: 
Koekemoer eraLl{T999t; 5:|Holtzman et al.|(T996); 6: [Donahue et al.|j2000f ; 8: |Donahue eral.. l ,2QQ7j ; 9; Edge et al., ( ,2010aj ; 10; .Edge et al.|p010b) ; 11: 
Sarazin et al.H1995) ; 12; [Taylor et al.nT999) ; 13; |Clarke et al.lpOOS) 



Observatory 


Instrument 


Mode/Config 


Band/Filter/Line 


Exp. Time 


Prog. / Obs. ID 


Obs. Date 


Program PI 


Reference 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


X-RAY Observations 


Chandra 


ACIS-S 


FAINT 


X-ray (0.5-7 keV) 


39.8 ks 


922 


2000 Jul 28 


McNamara 


1 




ACIS-S 


VFAINT 


X-ray (0.5-7 keV) 


52.8 ks 


6934 


2006 May 1 


Clarke 


New 




ACIS-S 


VFAINT 


X-ray (0.5-7 keV) 


60.9 ks 


7329 


2006 May 4 


Clarke 


New 


Far Ultraviolet / Optical / Near-Infrared Observations 


HST 


ACS 


SBC 


F150LP (FUV) 


8141 s 


11131 


2008 Jul 21 


Jaffe 


2 




STIS 


FUV MAMA 


F25SRF2 (Lya) 


1000 s 


8107 


2000 Jul 27 


O'Dea 


3 




WFPC2 




F410M 


2200 s 


6717 


1996 Jul 27 


O'Dea 


4 




WFPC2 




F450W 


2500 s 


6228 


1995 May 07 


Trauger 


5 




WFPC2 




F702W (R-band) 


2100 s 


6228 


1995 May 07 


Trauger 


5 




NICMOS 


NIC2 


F212N 


12032 s 


7457 


1997 Oct 19 


Donahue 


6 




NICMOS 


NIC2 


F160W (//-band) 


384 s 


7457 


1997 Dec 03 


Donahue 


6 


ESO VLT 


SINFONI 




/f-band 


600 s 


075.A-0074 


2005 Jul / Aug 


Jaffe 


7 


Mid- / Far-Infrared Observations 


Spitzer 


IRAC 


Mapping 


3.6, 4.5, 5.8, 8 /xm 


3600 s (each) 


3506 


2005 Nov 24 


Sparks 


8 




MIPS 




24, 70, 160 Atm 


2160 s (each) 


3506 


2005 Jun 18 


Sparks 


8 


Herschel 


PACS 


Photometry 


70, 100. 160 urn 


722 s (each) 


13421871(18-20) 


2009 Nov 20 


Edge 


9 




SPIRE 


Photometry 


250, 350, 500 ^m 


3336 s 


1342187329 


2009 Nov 30 


Edge 


9 




PACS 


Spectroscopy 


[O I] A63.18 fim 


6902 s 


1342187124 


2009 Nov 20 


Edge 


10 








[0 III] A88.36 Atm 


7890 s 


1342188703 


2009 Dec 30 


Edge 


10 








[N II] A121.9Atm 


7384 s 


1342188942 


2010 Jan 04 


Edge 


10 








[0 IB] XU5.52ij,m 


7382 s 


1342188704 


2009 Dec 30 


Edge 


10 








[C I] A157.74 ^tm 


6227 s 


1342187125 


2009 Nov 20 


Edge 


10 








[Si I] A68.47 fim 


11834 s 


1342210651 


2010 Dec 01 


Edge 


New 


Radio Observations 


NRAO VLA 




A array 


8.44 GHz 


15 min 


AR279 


1992 Nov 30 


Roettiger 


11 






A array 


4.99 GHz 


95 min 


BT024 


1996 Dec 7 


Taylor 


12,13 






A array 


1.3 GHz 


323 min 


BT024 


1996 Dec 7 


Taylor 


12,13 






A array 


330 MHz 


180 min 


AC647 


2003 Aug 18 


Clarke 


13 






B array 


330 MHz 


138 min 


AC647 


2002 Jun 10 


Clarke 


13 



|2009| l. New evidence suggests tliat star formation occurs when mul- 
tiphase clouds and filaments precipitate out of the ICM as its ce ntral 



entropy drops below a critical threshold (<; 20-30 keV cm , Voit 
et al II2008I ICavagnolo et al.||2008[ [Rafferty et al.||2"008l|Sharma 



et aL|2011[[Gaspari et al.|20Tl[ i. 



But while the pathway of entropy loss from hot ambient 
medium to cold star forming clouds may be strongly influenced by 
AGN feedback, it is not known whether residual cooling persists 
at constant low levels, or in elevated episodes anti-correlated to the 
AGN duty cycle (e.g., [O'Dea et al.|2010[[Tremblay|20TT] l. More- 
over, models invoking radio-mode feedback to quench cooling flow 
signatures such as star formation must be reconciled with evidence 
that, in several systems (e.g., Abell 1795 and Abell 2597, |0'Dea| 
|et al.|[2004l >, the propagating radio source may trigger compact, 
short-duration starbursts as it propagates amid a dense medium 



(e.g.,IElmegreen & Elmegreen|1978, Voit|1988 


iDe Young 1989 


McNamara & O'Connell 1993|fo'Dea et al.|200^ 


||Batcheldor et al. 


2007;iHolt et al.|2008,|20l l|l. The picture is further complicated by 



the role played by (e.g.) thermal conduction and cold accretion sce- 



narios like gas-rich mergers, whose importance relative to residual 
hot-mode ICM cooling has driven a debate that is still unsettled af- 
ter three decades (e.g. [Holtzman et al.||1996[ [Sparks et ar]|1989| 
[2012| >. Advances in understanding largely rely upon multiwave- 
length data that sample all temperature phases of the ISM in CC 
BCGs, the transport processes between these phases, and their as- 
sociated mass and energy budgets. 



1.2 The brightest cluster galaxy in Abell 2597 

To that end, in this paper we present a multiwavelength study 
of the archetypal cool core cluster Abell 2597. New Chandra X- 
ray and Herschel FIR observations, combined with a vast suite of 
archival data, enable a radio-through-X-ray study of the ICM/ISM 
on the scale of its AGN feedback interaction region (i.e., its central 
<; 30 kpc X-ray cavity network). Abell 2597 (hereafter A2597) is 
an Abell jAbell[[T958l [Abell et al.[[T989l l richness class galaxy 
cluster with an X-ray surface brightness profile that is sharply 
peaked about its centrally dominant elliptical BCG at redshift 
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z = 0.0821 ( jVoit & Donahue||l997^ . The galaxy is host to one of 
the nearest known compact steep spectrum (CSS, [O^Dea|1998| ra- 
dio sources, PKS 2322-122, which exhibits a compact (10 kpc) and 
bent Fanaroff-Riley class 1 (FR 1) morphology at 8.4 GHz ( Fana roff] 
& Rileyl 1974) [Wright & Otrupcek| 1990, Griffith & Wright 1994| 



O'Dea et al.|1994bl|Sarazin et al.|1995> . 

A2597 is one of the only CC clusters with a convincing FUSE 
detection of [O VI]A1032 A emission stemming from a 10^ - lO' 
K gas component ( Oegerle et al.'2001 ), and high resolution XMM- 
Newton X-ray specroscopy reveals weak Fe XVII features amid a 
soft X-ray excess ^Morris & Fabian|2005^ . It is therefore one of the 
best known candidates for harboring a moderately powerful resid- 
ual cooling flow with a mass deposition rate of 90 ± 15 Mq yr"' 
within a 100 kpc cluster-centric radius (iMorris & Fabian 2005|. 
The BCG harbors a 1.8 ± 0.3 x l O'' cold H2 component, in- 



Salome & Combes 



ferred from CO observationQ(e.g., Edge|20oT 
|2003 |), as well as a young stellar component cospatial with an Ha 



bright filamentary emission line nebula (e.g. Heckman et al. 



[Wr& Donahue 1997; Koeke moer etal.|1 999; O'D ea et al. 



1989 



2004 



|Oonk et al._2011j . These features reside amid a network of promi 
nent X-ray cavities ( [McNamara et al.|[200H [Clarke et ar]|2005| l, 
making A2597 an ideal subject for studies of AGN/ISM interac- 
tions ( [Tremblay et al.|2012a[ l. For these reasons and others, A2597 
enjoys a long history of cross-spectrum analysis in the literature 
(X-ray —[Crawford et al.'1989'.'Sara zin et al.|T995|[Sarazin & Mc- 
Namara] 1997, McNamara et al. 2001 [ [Clarke et al.|2005[ [Morris 
& Fabian 2005 , UV/optical — McNamara & O'Connell 1993;'De 
Young 1995 , Voit & Donahue|1997 , Cardiel et al. 1998 , McNamara 



etal.|1999 



Koekemoer et al.|1999 



jOegerle et aL|2001 



_2004MJaffe et al.«2005','Oonk et al.2010, 20irjlR 



O'Dea et alT] 
]McNamara' 



& O'Connell 1993 , Voit & Donahue 1997; Mc Namara et al. 1999 
koekemoer et al. 1999, O'Dea et al . 2004, Jaffe'eFair 2005 , Don- 
ahue et al..2007 201 1 , sub-mm — Edge 2001 . Salome & Combes 



2003^ 



Edge et al.[[2010b|a[ radio —[O'Dea et al.[[1994a[|Sarazin 



erar][1995| [Taylor et al.][1999l|Pollack et al.[|2005| [Clarke et al 



2005b 



In Section 2 we describe the new and archival multiwave- 
length observations used in this paper. In Section 3 we use the 
new Chandra X-ray observations to assemble an AGN feedback 
energy budget by using the kpc-scale cavity network as a lower- 
limit calorimeter to the AGN kinetic energy input. In Section 4 we 
present new Herschel FIR data and use it to place constraints on the 
residual cooling flow. These results are used to frame a discussion 
on star formation, which we present in Section 5. A concluding dis- 
cussion and summary are provided in Sections 6 and 7. Through- 
out this work, we adopt Ho = 71 h^l km s"' Mpc"', Q,m = 0.27, and 
Qa = 0.73. In this cosmology, 1" corresponds to ~ 1.5 kpc at the 
redshift of the A2597 BCG (z = 0.0821). This redshift corresponds 
to an angular size distance of Da ~ 315 Mpc and a luminosity dis- 
tance of Dl ~ 369 Mpc. 



2 OBSERVATIONS & DATA REDUCTION 

Table [T| contains a summary of all new and archival observations 
used in this paper. We refer the reader to the publications listed 
in column (9) for technical information pertaining to the archival 
observations. The new Chandra X-ray observations, totalling 150 

^ Updated CO-infeired molecular hydrogen mass is from P. Salome (2012, 
pvt. comm.), using unpublished IRAM observations. We discuss this in Sec- 
tion 4. 



ksec (128 ksec flare-free) in combined effective exposure time, 
were reduced using the standard CIAO v4.2 threads with v4.3.1 
of the calibration database. More details can be found in |Tremblay[ 
[et al.[ ( [2()12 a), which presents the spatial and spectral analysis of 
these new observations. 

Herschel Space Observatory observations of A2597 were ob- 
tained in November 2009 with the Photodetector Array Camera 
& Spectrometer (PACS), as well as the Spectral and Photomet- 
ric Imaging Receiver (SPIRE). These observations were part of an 
Open Time Key Project (OTKP) investigating the FIR line and con- 
tinuum properties of a sample of 1 1 BCGs in well-known X-ray 
and optical line selected cool core clusters (PI: A. Edge, 140 hrs). 
Preliminary results for A2597 have been published in the works 
submitted as part of the science demonstration phase ( [Edge et aTj 
[2010am . The data were processed with the Herschel Interactive 
Processing Environment (HIPE) package version 7.1.0. See [Edge] 
[etaLl ( [2010a|b) , [Mittal et al.[pOT ri, and Oonk et al. (2012, in prep) 
for details on the data reduction. The present paper uses PACS 
and SPIRE data to fill in critical gaps in the FIR SED for A2597, 
enabling constraints on dust masses and temperatures. The [O I] 
A68.4 fim PACS observation, because it is at high spectroscopic 
resolution, is also used to place constraints on the cold gas kine- 
matics. 



3 AGN FEEDBACK ENERGY BUDGET 

Kiloparsec scale X-ray cavities can be used as lower-limit 
calorimeters on the kinetic energy associated with AGN outflows, 
and the duty cycle over which these outbursts occur (e.g.,'Churazov 



etal. 



2006 



2002 ; Blrzan et al. 2004 , Dunn & Fabian 2006 ; Rafferty e taL 
McNamara & Nulsen||2007| |2012( ^Fabian 20 12||l. A study 



of the original 40 ksec Chandra observation by M cNamara et al.[ 
pOOT) and [Clarke eTaL] pOOSt (hereafter MOl aiid C05, respec- 
tively) revealed western and northeastern ghost cavities in A2597. 
The new 150 ksec Chandra data reveals the cavity network to be 
more extensive than originally thought, changing the interpretation 
of the AGN feedback energy budget and outburst history. We dis- 
cuss these two new results in this section. 



3.1 The X-ray cavity network 

In Fig.[T]we present the combined 150 ksec 0.5-7 keV X-ray data 
from Tremb lay et al.[p012a^ . All panels are aligned with an iden- 
tical 50" X 50" (75 kpc x75 kpc) field of view (FOV) centered on 
the X-ray centroid at RA=23h 25m 19.75s, Dec =-12° 07' 26.9" 
(12000). East is left, north is up. The top left panel shows the ex- 
posure corrected, merged data adaptively smoothed with a variable 
width gaussian kernel whose radius self-adjusts to match the lo- 
cal event density. Black surface brightness contours are overlaid 
to make individual features easier to view. The innermost contour 
marks a flux of 4.6 x 10"' photons sec"' cm"^ pixel"', and the con- 
tours move outward with a decrement of 2.0 x 10"** photons sec~' 
cm"^ pixel"' . 

In the top right panel of Fig. [T] we overlay the adap- 
tively smoothed X-ray contours on the Hubble Space Telescope 
(HST) WFPC2 /;-band observation of the A2597 BCG (shown in 
red/orange). Note that the anisotropy of the X-ray emission is con- 
fined to the scale of the BCG, while the outermost regions assume 
a smoother, more elliptical shape. The major axes of the X-ray and 
BCG stellar isophotes are aligned. 

In the bottom two panels of Fig. [TJwe show the same 0.5-7 
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Figure 1. (top left) Exposure-conected 0.5-7 keV image of the three merged Chandra observations (see Table 1). The data have been smoothed with an 
adaptive gaussian kernel. Black contours are overlaid to better show the spatially anisotropic nature of the emission. The innermost contour marks a flux of 
4.6 X 10"' photons sec"' cm"^ pixel"' , and the contours move outward with a flux decrement of 2.0 X 10"** photons sec"' cm"- pixel"' . (top right) The same 
0.5-7 keV X-ray contours (in green), overlaid on the broad-band //577WFPC2 F702W (S-band) exposure of the Abell 2597 brightest cluster galaxy (shown in 
orange). Note that the anisotropy of the X-ray emission is largely confined to the inner regions of the galaxy, (bottom left) Unsharp mask of the X-ray data, 
made by subtraction of a 20" Gaussian smoothed version from the adaptively smoothed version. 330 MHz VLA contours have been overlaid in green, while 
the 1.3 GHz radio contours are plotted in blue, and the 8.4 GHz radio contours appear in black, (bottom right) Unsharp mask made by subtracting the same 20" 
Gaussian smoothed version of the image from a 5" Gaussian smoothed version, then dividing by the sum of the two images. The major (lOcr excess or deficit) 
morphological features which will be the subject of further spatial analysis in this paper have been labeled 1-6. The color schemes of the two bottom panels 
show regions of X-ray surface biightness excess in red/orange, while deficits (cavities) appear in black. All panels share an identical field of view, centered at 
RA=23h 25m 19.75s, Dec =-12° 07' 26.9" (J2000). 
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keV data processed in two ways: (1) a highly processed unsharp 
mask made by subtracting a 30" gaussian smoothed image from 
the adaptively smoothed image shown in Fig. [TJi, and (2) a more 
conventional unsharp mask made by subtracting a 20" gaussian 
smoothed image from a 5" (non-adaptive) gaussian smoothed ver- 
sion. The subtracted data is then divided by the sum of the two 
images. Regions of X-ray surface brightness excess over the sub- 
tracted smooth background appear in white, while deficits (cavities) 
appear in black. 8.4 GHz, 1.3 GHz, and 330 MHz Very Large Ar- 
ray (VLA) radio contours from | Sarazin et al^ ( |1995[ l, [Taylor et al.| 
([l|99), and [Clarke et al] \20Q5) are overlaid on the bottom left 
panel. See section 3.2 of |Tremblay et al.] ( |2012a| l for a discussion 
of the strong radio/X-ray spatial correlations that are evident in this 
figure. We note that these unsharp mask edge enhancement meth- 
ods (particularly method 1) are inherently noisy, and introduce ar- 
tifacts that complicate quantitative morphological analysis. Signif- 
icance of the observed features must therefore be estimated from 
the un-processed data. We show these mostly as viewing aids for 
the X-ray cavities discussed in this section. We do note that the 
processed images shown here do not significantly differ in appar- 
ent morphology from the un-processed raw Chandra data that can 
be seen in Fig. 2 of |Tremblay et al.| ( |2012a i 

In the bottom right panel of Fig.fllwe label the six features that 
are associated with ^ lOa deficits or excesses relative to the local 
mean. The significance of these features was estimated by com- 
paring unsmoothed counts in like sized regions at the same cluster- 
centric radius. Throughout this paper we refer to these features with 
a common label and name. These are: (1) the "MOl western ghost 
cavity", (2) the "COS X-ray tunnel", (3) the "MOl northern ghost 
cavity", (4) the "eastern ghost cavity", (5) the "cold filament", and 
(6) the "COS filament base cavity". For clarity, these names and la- 
bels will be used consistently throughout this paper. Feature (2), the 
"cold filament", is discussed in Section 5 of |Tremblay et al.| ( |2012a^ 
in the context of AGN/ISM interactions and multiphase gas dredge- 
up by the radio source. 

While (for "historical" reasons) we label features (1) and (2) 
independently, our deeper data make it clear that western ghost cav- 
ity described by MOl and COS is part of a larger "teardrop" shaped 
cavity ~ 2S kpc in projected length (COS originally suggested this 
in their discussion of an "X-ray tunnel"). This changes the inter- 
pretation of the AGN outburst history (relative to the conclusions 
drawn in MOl), which we will discuss below. So as to enable com- 
parison with past papers (e.g., M01; [Birzan et al.|2004| [Rafferty| 
|et al. [2006) 1 that have treated the MOl cavity as a separate feature, 
we maintain independent labels (1 and 2) for the cavity and tunnel. 
We stress that these features are part of one larger cavity, which we 
will label "1+2" and call the "western large cavity" for the remain- 
der of this paper. 



3.2 Age dating the X-ray cavities 



We adopt the simple model used by |Birzan et al.|p004[ > and |Raf-| 
[ferty et aT] ( |2006| l in estimating rough ages for the X-ray cavi- 
ties, assuming their density is very low relative to the ambient gas 
density. Three timescales can be considered, the most simple of 
which assumes that the bubb le rises in the plane of the sky at the 
sound speed c, ~ kT/nntp, where = 0.62 is the mean molecular 
weight in units of the proton mass, nip. If the bubble is at projected 
cluster-centric radius R, then this sonic rise timescale is given by 



As the initial stages of cavity inflation are thought to be supersonic, 
followed by subsonic buoyant rise, this simple approach may best 
reflect an average of the two phases. Alternatively, if the supersonic 
inflation period is a negligible fraction of the cavity's age and drag 
forces limit the cavity's terminal velocity v,, the buoyant timescale 
can be written as 



R 

V, 



■R 



AC 
2gV' 



(2) 



where g is the local gravitational acceleration, V is the cavity vol- 
ume, A is its cross-sectional area, and C = 0.7S is the drag coef- 
ficient adopted from' Churazov et^r] ( [200l| by (e.g.) |Birzan et al.] 
(20041. Finally, cavity ages can be constrained by estimating the 
time required to refill the volume displaced by a bubble of radius r 
as it rises a height equal to its diameter. 



'refill — 2R 



GM«R) 



(3) 



where M(< R) is the mass enclosed within a sphere of radius R. 

[Birzan et al.| ( |2004) made all three estimates for the NE and 
W ghost cavities in A2S97 (features 1 and 3 in Fig.[T] bottom right 
panel), using the early short (40 ksec) Chandra observation. So as 
to provide an independent check of their results using the deeper X- 
ray data, we follow their procedure almost exactly, and repeat their 
calculations for the MOl ghost cavities, the western large cavity 
(features I and 2), as well as the newly detected cavities (4) and (6) 
(as labeled in Fig.[T](. Like |Birzan et al. ( 2004 1, we adopt the A2S97 
BCG stellar velocity dispersion ofa~ 224 ± 19 km s"' from Smith 



|et al.|p990) in our calculation of the local gravitational acceleration 



2a^ 
R 



(4) 



assuming that the galaxy is an isothermal sphere. In |Tremblay et al.| 
( |2012a| l we calculate a local value for g by estimating the total en- 
closed gravitating mass from a /3-model fit to the X-ray surface pro- 
file. For this discussion however, it is sufficient to adopt the Bir zan] 
|et al.| ( (2004) method for estimating g (the two methods turn out to 
be roughly consistent at the radius of the X-ray cavity network any- 
way). Unlike [Birzan et ari ( |2004| , we use the kT inferred from the 
cavity positions on the X-ray temperature map presented in |Trem-[ 
[blay et ari ( [20I2al l to calculate the sound speed in the X-ray gas. 
We use the 2D temperature map in lieu of the ID radial temperature 
profile because A2S97 is azimuthally anisotropic in X-ray temper- 
ature (not to mention surface brightness) on these scales. Further- 
more, in calculating the pV work associated with each cavity, we 
use the projected density profile to estimate the pressure nkT at the 
cluster-centric radius of each cavity. 

The results of these calculations are given in Table 2. Our 
findings are roughly consistent with those of 'B irzan et al.[p004[ ). 
For example, in age dating the MOl western ghost cavity we find 
tc, « 27 Myr, /buoy ~ 88 Myr, and /refill ~ 66 Myr, while [Birzan et al.[ 
( [2004) 1 finds 26, 66, and 86 Myr, respectively. None of these esti- 
mates account for projection effects, and all assume that the bubble 
rises purely in the plane of the sky. This results in underestimation 
of cavity age by generally less than a factor of two (Tjirzan et aT] 
[2004|[Rafferty et al.|2006l[McNamara & Nulsen|2007y 



tc, ~ R/c, 



(1) 



3.3 Cavity heating energy reservoir 

In column (9) of Table 2 we calculate the mean instantaneous power 
of each cavity Pcav, 
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Table 2. Spatial properties and energetics of the X-ray cavities. (1) Label of morphological feature that corresponds to that assigned in the 
lower right panel of Fig.[T] (2) name given to the corresponding feature; (3) projected radial distance (length) of the feature from the radio 
core to the estimated center (edge) of the feature; (4) estimated radius of the feature; (5) estimated work associated with cavity assuming 
subsonic inflation; (6) age of the cavity if it rises at the local sound speed; (7) buoyant, subsonic cavity rise time; (8) time needed to refill 
the displaced cavity volume; (9) X-ray cavity power, assuming the cavity is filled with relativistic plasma. See Section 3 for more details on 
these calculations. 







R(D) 


r 


pV 




^buoy 


'refill 


^cav 


Label 


Name 


(kpc) 


(kpc) 


(xlO^^ ergs) 


(xlO' yr) 


(xlO^ yr) 


(xlO' yr) 


(Xl0**2ergs-') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


1+2 


Western "large cavity" 


9 


9 


35.9 


1.0 


1.4 


5.6 


170.6 


3 


MOl northern "ghost" cavity 


21 


6.6 


7.0 


2.7 


6.1 


7.3 


16.5 


4 


Eastern "ghost" cavity 


35 


3.6 


0.79 


3.8 


17.8 


6.9 


1.05 


6 


Filament base cavity 


9 


2.3 


0.30 


1.1 


2.7 


2.8 


1.73 



4pV 



(5) 



where (t) is the average of the cavity ages listed in columns 6, 7, 
and 8. We have assumed the cavities are filled with a relativistic 
plasma, such that their enthalpy can be approximated as 4pV (e.g., 
[McNamara & Nulsen|2007l l. 

The total sum of all cavity thermal energies (pV) listed in col- 
umn (5) is 4.4 X 10^* ergs. This serves as a rough, lower-limit es- 
timate on the kinetic energy injected by the AGN into the ambient 
X-ray gas during the (a) current AGN episode or perhaps (b) dur- 
ing the past two or more episodes of activity. This is close to the 
inferr ed mechanical energ y of the central 8.4 GHz radio source, 
which 



Sarazin et al. 



1995 I estimated to be 9 x 10^ ergs. In princi- 



ple, the roughly estimated mechanical input from the radio source 
could be capable of accounting for the energy budget inferred from 
the X-ray cavities, ignoring timescale arguments. Whether or not 
these cavities have been produced by one or more AGN episodes is 
the subject of the following section. 



3.4 The AGN duty cycle and heating timescale budget 

A network of multiple cavities found at varying cluster-centric radii 
(as we find in A2597) may be produced with an episodically vary- 
ing AGN, transitioning between either "on" and "off" or high and 
low modes. Cycling times between the triggering of radio activity, 
the onset of quiescence, and the subsequent re-ignition of activity 
are typically estimated to be on the order of 10^ - 10* yr in general 



(e.g.,|Parma et al.|1999[[Best et al.|2005||Shabala et al.|2008||TYen>l 



|blay et al.||2010 i, and synchrotron losses limit radio source life- 
times to ~ 10 yr, unless there has been re-acceleration of the elec- 
tron population. Alternatively, a steady-state AGN with a duty cy- 
cle near 100% can also produce a series of discrete cavities (rather 
like a fish tank aerator or a dripping faucet, e.g., [Peterson & Fabian| 
|2006[ [McNamara & Nulsen 2007 1. In this section we attempt to 
distinguish between these two possibilities. 

Based on synchrotron loss timescales, [Clarke et al.|(2005[ l es- 
timated the minimum energy, lower-limit lifetime of the 330 MHz 
source ( green contours in the bottom left panel of Fig. [TJ to be 
in the range of 1^30 mhz ^ 8 x 10* yr (if the spectral index is steep 
down to 10 MHz) to ;33o mhz ^ 5 x 10'' yr (if the spectral index flat- 
tens beyond 330 MHz). This model assumes equal cosmic ray and 
magnetic field energy densities (equipartition), a steep spectral in- 
dex of -2.7 between 1.3 GHz and 330 MHz, and that the source is 
a uniform prolate cylinder with a filling factor of unity. This model 
yields a minimum energy magnetic field strength of 29 pG and a 



non-thermal pressure of 5 x 10"" dyn cm"^. We note that equipar- 
tition may be a less-than-ideal assumption in this case (e.g., Fabi^ 
[et al.|2 002 ). Furthermore, deeper multiband radio observations are 
required to better assess the spectral index of the source across the 
western large cavity. 

The lower end of the 330 MHz lifetime range quoted above 
(fjsoMHz ^ 8 X 10* yr) is significantly shorter than the estimated 
buoyant rise time of the western large cavity it fills (see Table 2, 
column 7). If we assume that this is the actual age of the source, 
and further assume that the estimated age of the western large cav- 
ity is roughly correct, then we require (a) in situ re-acceleration 
of electrons in the cavity or (b) a new episode of activity to feed 
the plasma into an already pre-established cavity (which, if empty, 
should collapse on a sound crossing time <; lO' yr). As we will 
elaborate upon below, there is little supporting evidence favoring 
either of these scenarios. 

A "radio source younger than the cavity" scenario could also 
be possible if adiabatic and/or inverse Compton losses dominate 
over the synchrotron processes, which would result in the radio 
source fading on a timescale shorter than the synchrotron lifetime. 
The importance of these processes relative to synchrotron losses de- 
pends on the B field strength and unknown details of how the radio 
source fills the bubble. Inverse Compton losses become more effi- 
cient than synchrotron losses when the B field is extremely weak 
(S < 1 /iG, e.g., [Fabian et al. [2002) 1, and adiabatic losses are not 
likely a limiting factor because the western large cavity is likely 
rising at less than the sound speed (note the absence of any de- 
tected fast shocks in the X-ray data, e.g., [Tremblay et al.|2012a] l. 
We therefore consider this scenario unlikely. 

The upper end of the minimum-energy synchrotron lifetime is 
f330 MHz ^ 5 X 10' yr (which itself is merely a lower limit). The 330 
MHz radio source could therefore easily be roughly the same age as 
the cavity it fills. This is the simplest scenario, as it is assumed that 
the propagation of the radio source amid the X-ray gas is the mech- 
anism that excavates the cavity. Even if the equipartition and spec- 
tral index assumptions mean that this estimated lifetime is wrong, 
pressure equilibrium between the radio source and the ambient hot 
gas can be assumed to estimate a very similar ^ lO' yr minimum 
age (e.g., |Fabian et al. | i2002j l. Considering the high uncertainties 
associated with age dating both cavities and radio sources, (not to 
mention projection effects), we suggest that the western large cav- 
ity was created by the current, ongoing episode of AGN activity. 

Considering the above, we cannot rule out the possibility that 
the AGN in A2597 is on nearly ~I00% of the time. The oldest 
cavity (i.e., feature 4, with an estimated age of 1.8 x 10** yr) may 
indeed be associated with a previous epoch of AGN activity, but it 
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could also be a detached bubble excavated by the current episode 
(the available data does not permit us to discriminate between these 
possibilities). We note that the position angle of the 8.4 GHz ra- 
dio lobes are significantly offset from the western large cavity, 330 
MHz and 1.3 GHz radio source (Fig. [T), an d the VLBA small- 
scale (~ 50 pc) jet axis ( [Taylor et al.||1999) . We find it unlikely 
that the offset is due to the 8.4 GHz source being associated with 
a previous epoch of activity, considering its steep spectrum ( |0'Dea| 
|et al.||1994a| [Clarke et al. 20051, estimated minimum-energy age 
of > 5 X 10 yr (Sarazin et al.,1995^ and alignment of the VLBA 
(< 50 pc) and 330 MHz (> 25 kpc) major axes ( [Tremblay et al.| 
[2012a| l. Several past studies have suggested that the bend is due 
to interaction of the radio source with the ambient dense gaseous 
medium, either by deflection or back-flow along pre-established 



pressure gradients (e.g.,[Sarazin et al. 


1995 


Koekemoer et al.|1999| 


[O'Dea et al.|2004| [Clarke et al.|2005 


lOon 


k et al.|2010|l. 



the corresponding range of cooling luminosity is roughly Lcooi 
(1 - 4) X lO*"* ergs s ~'. T he rough sum of mean instantaneous cavity 
powers (see Section : 



3.3 



Jand column 9 of Table 2) is P,.av -2x10** 
ergs s"' , suggesting that, in principle, there is enough enthalpy asso- 
ciated with the current cavity network to offset the bulk of radiative 
losses and effectively quench the cooling flow. 

The same conclusion holds if we estimate the predicted cool- 
ing flow luminosity directly from the new Chandra data, us- 
ing lower-limit temperatures derived from spectral fits (i.e., the 
MKCFLOW model A:?]™ values in Table 1 of [Tremblay et al.[2012al l. 
This is a rough, order-of-magnitude assumption which neglects 
many important considerations such as how this enthalpy might be 
dissipated and distributed in the ambient gas. Our finding that the 
X-ray cavities are associated with enough energy to inhibit the pre- 
dicted classical cooling flow luminosity is consistent with previous 
results jBIrzan et al.|2004|[Dun"n & Fabian|2006t . 



3.5 Black hole accretion rate 

If the AGN is indeed nearly steady state, the central black hole (BH) 
requires a stable supply of gas from the ambient accretion reservoir. 
Assuming a mass-to-energy conversion efficiency of e = 0. 1 ^Wise[ 
[et al.|2007[ l, and assuming that the energy associated with the X-ray 
cavities is provided by the AGN, the sum of mean instantaneous 
cavity powers (Pc„) requires a time averaged BH accretion rate of 



'0.003-0.03 



OA 



Mq yr 



(6) 



where c is the speed of light. Using the Mbh-ct relation ( Magorrianj 
[et al.|1998[[Ferrarese & Merritt|2000[ [Gebhardt et al.|2000| l with 
the ^-band host luminosity and stellar velocity dispersion yields a 
rough BH mass estimate of ~ 3 x 10* Mq, for which the c orre- 
sponding Eddington accretion rate would be ~ 10 Mq yr"' iRaf- 
[ferty et al.|2006^ . This suggests that the accretion rate, if steady, is 
strongly sub-Eddington, consistent with past results even for very 
powerful AGN (e.g.. Hydra A, [Wise et al.[[2007l >. If the AGN is 
instead strongly variable or episodic, then gas transport to the BH 
could be non-steady and the actual accretion rate could vary. 



3.6 There is enough energy to quench a classical cooling flow 

We calculate the classical X-ray derived cooling time for A2597, 
5 nkT 



(7) 



where n and kT are the gas density and temperature profiles ob- 
tained by [Tremblay et al.] ( [2012a^ , and A is the T > 0.02 keV por- 
tion of the cooling function from'Sutherland & Dopita[ ( [1993[ l, us- 
ing the generalization from Tozzi & Norman 200 1[ 

We find that the cooling time at the <; 30 kpc outermost radius 
of the X-ray cavity network is ~ 300 Myr. The instantaneous cool- 
ing luminosity associated with the uninhibited cooling flow (with 
no heating) can be calculated from the classical uninhibited mass 
deposition rate Mcooi using 



5Mc, 



-kT 



(8) 



2 /iOTp 

where Tvi, is the virial temperature of the cluster. Classical mass 
deposition rates for A2597 range from ~ 100-500 Mq yr"' over a 
range of cluster-centric radius spanning ~ 30- 100 kpc (e.g., [Allen] 
[et al.|2001[ [McNamara et al.|2001[ [Morris & Fabian|2005[ >. Tak- 
ing a cluster virial temperature of ^Jvir = 3.5 keV, we find that 



4 NEW EVIDENCE SUPPORTING A RESIDUAL 
COOLING FLOW MODEL FOR A2597 

In the previous section we estimated that the X-ray cavity network 
is associated with enough enthalpy to (in principle) replenish a sig- 
nificant fraction of the radiative losses associated with a classical 
cooling flow. However, the warm and cold gas phases in the A2597 
BCG strongly suggest that some cooling from the ambient X-ray 
atmosphere has managed to persist, even amid AGN heating. We 
discuss this possibility here. 



4.1 How might cooling persist amid AGN feedback? 

As mentioned in Section 1, results from XMM-Newton and FUSE 
are consistent with a moderately strong residual (sometimes called 
"reduced") cooling flow with a mass deposition rate of ~ 30± 15 
Mq yr"'at ~ 30 kpc to 90 ± 15 Mq yr"' at 



100 kpc ( [Oegerle 
[et al.|200T| [Morris & Fabian 2005V The cooling luminosity asso 
ciated with this residual cooling flow is Lcooi.Rcsid ~ (0.3 - 1) x lO'*'* 
ergs sec"'. As with 15-30% of CC BCGs (depending on the sam- 
ple, e.g., |Salome et al.[2006| >, the A2597 BCG harbors a substantial 
1.8 ±0.3 X 10^ Mq reservoir of cold H2 (inferred from CO obser- 
vations) within its central 30 kpc (Edge 2 0()T] [Salome & Combes[ 
[2003[ P. Salome private communication, 2012). If the bulk of the 
mass budget for this cold gas is supplied by the residual cooling 
flow (rather than e.g., a merger), then AGN feedback cannot estab- 
lish an impassable entropy floor, and some cooling to < 100 K must 
be permitted even if there is enough energy in principle to quench 
the classical cooling flow. We consider four possibilities: 

(i) Low levels of residual cooling persist even while AGN feed- 
back is heating the ambient environment; 

(ii) cooling occurs in a spatially structured manner, away from 
those regions which are being locally heated by (e.g.) buoyant cav- 
ities and sound waves; 

(iii) cooling occurs in episodes which correspond to the times 
when the BH is inactive, or transitioning between active and inac- 
tive states or between high and low modes; 

(iv) the substantial cold gas reservoir in the A2597 BCG stems 
from something other than a cooling flow, i.e. one or more gas rich 
mergers. 

These are not the only possibilities, and several or all of the four 
scenarios listed above could work in tandem. 
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As discussed in Section [14| we cannot rule out the possibility 
that the AGN duty cycle in A2597 is close to 100%. Even if this is 
not the case, and the AGN is episodic over lO' yr periods, a cav- 
ity's enthalpy reservoir should be dissipated as heat in the ISM/ICM 
on a timescale of order the cavity lifetime ( Churazov et al.|2002| 
[Reynolds et~ar][2002i [Birzan et alT][2004,i Dunn & Fabi an||2006X 
The net effect of this time buffer could be that ISM/ICM heating 
rates are roughly constant in time for cavity heating models, even if 
the AGN itself is episodic over these timescales. If we assume that 
this is correct, and further assume that most cold baryons in the 
nucleus cooled from the hot atmosphere, then we might conclude 
that low levels of residual cooling must persist at a roughly con- 
stant rate. This would be consistent with scenario (i) or (ii) above. 
In this case, a near-constant cooling rate would funnel a steady flow 
of gas to the BH accretion reservoir, fueling a long-lived rather than 
episodic AGN. This would be consistent with our suggestion that 
the AGN in A2597 is approximately steady-state. We would require 
inhomogeneous heating and cooling (e.g., scenario ii, iii and/or iv) 
to allow cooling to re-trigger the radio source. 

The notion that cooling could proceed in a spatially discrete 
and structured manner (scenario ii) is not controversial, and was 
predicted even by the earliest cooling flow models (see e.g., the re- 
view by |Fabiari|1994 ). While long-lived sound waves could theoret- 
ically heat the outermost regions of cool cores in a more homoge- 
nous manner, heating by short-lived X-ray cavity enthalpy dissipa- 
tion is expected to be spatially confined to the scale and location 
of the cavity, regardless of the Reynolds number of the plasma (see 
e.g., the review by McNamara & Nulsen 2007). In this case, one 
would expect cooling to proceed in the regions unaffected by local 
heating. Recently, deep X-ray observations of Perseus and A2146 
by (respectively) [Fabian et al.| ( |2011 t and Russell et al. (2012) re- 
vealed bright X-ray filaments which alone could be associated with 
a significant fraction (< 80%) of the total inferred cooling rate. In 
Section 5 of this paper, we will show a similar result for A2597. 

On the other hand, [O'Dea et al.[ ( |2010| l studied a sample of 
7 CC BCGs selected on the basis of an IR excess associated with 
elevated star formation rates, and found that each tended to possess 
a weak, compact radio source. The combination of higher SFR and 
lower radio power may be consistent with a scenario wherein a low 
state of AGN feedback allows for increased residual condensation 
from the ambient X-ray atmosphere, accounting for the elevated 
star formation rates. This would be consistent with scenario (iii), 
though A2597 does not appear to be a strong candidate. We stress 
that scenarios (i), (ii), and (iii) are not contradictory or competitive 
with one another, and each could play a simultaneously important 
role. The same is true for possibility (iv), which we investigate in 
more detail below. 



4.2 New results on cold gas kinematics — discriminating 
between iiot and cold accretion scenarios 

The new Chandra and Herschel data can be used to better under- 
stand the origin of the 1 .8 ± 0.3 x lO' Mq cold reservoir in A2597. 
Relative to BCGs in non-cool core clusters, CC BCGs are far more 
likely to harbor central emission line nebulae, cold gas components, 
and radio sources (e.g., [Mittal et al.|2009[[Hudson et al.|20T0| and 
references therein). This strong circumstantial evidence suggests a 
causal connection between these warm/cold phenomena and resid- 
ual condensation from the ICM, but the literature has so far not 
reached consensus as to the relative importance of gas-rich major 
and minor mergers in this relationship (e.g., I Sparks et al.[[2012| ). 
[Koekemoer et al.[ jl999[ l, for example, suggest that some physical 
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[km/s] 

Figure 2. A comparison of the Herschel [O I] A63/jm line profile, in black, 
with the TRAM 30 m CO(2-l) rotational line profile in blue. The [O I] line is 
centered at the systemic BCG redshift, and has an intrinsic FWHM of 405 it 
55 km sec"' with a possible asymmetric red velocity excess offset from the 
systemic velocity by ~ +250 km sec"'. The CO(2-l) FWHM is slightly 
narrower (300 it 50 km sec"' ), and may also exhibit a slight red asymmetry. 
The spectral resolution of the CO data is 45 km s"' . Both spectra have been 
normalized to their maximum to enable comparison of the line profiles. 
Intensities are therefore arbitral^. 



features of the A2597 emission line nebula are consistent with what 
might be expected from a gas-rich minor merger The BCG stellar 
isophotes do not show any obvious signs of disturbance, so a recent 
major merger is likely ruled out, at least within the past several dy- 
namical times. 

In Fig.|2]we plot the [O I] A63 pm line profile from the Her- 
schel PACS spectroscopy jEdge et al.|[20I0b) in black, over the 
C0(2-l) rotational line profile (in blue) from the IRAM 30 m ob- 
servations of A. Edge and P. Salome (pvt. comm. 2012, data pre- 
viously unpublished). [O I] A63/xm is one of the primary atomic 



coolant lines for T ^ 40 K gas i Kaufman et al. 1999 1, and C0(2-l) 
can be considered a tracer of the molecular hydrogen reservoir Fol- 
lowing Solomo n et al.[ jI99 7 1 and adopting a standard Milky Way 
H2 mass-to-CO luminosity ratio of a = 4.6 Mq (K km s"' pc^)"', 
the CO line luminosity of Lqq = 3.9 x 10^ K km s"' pc" converts to 
a CO-inferred molecular hydrogen mass of Mh, = 1 .8 ± 0.3 x lO' 

Mq. 

[O I] was detected at S/N^30 but unresolved (the Herschel 
beam at this wavelength is ~ 9" ~ 13.5 kpc). The [O I] line is cen- 
tered at the systemic BCG redshift, and has an intrinsic FWHM of 
405 ± 55 km sec"' with a possible asymmetric red velocity excess 
offset from the systemic velocity by ~ +250 km sec"' . The C0(2-l) 
FWHM is slightly narrower (300 ± 50 km sec"' ) though the profile 
may also exhibit a scaled down signature of the same red asymme- 
try ([Edge etal.|2010b|). The IRAM 30 m C0(2-l) beam is 13" ^ 20 



kpc, larger than the Herschel 63/im beam, so a red wing arising 
from gas near the center would be more diluted in the C0(2-l) pro- 
file. A very marginal {2.6a) CO(l-O) detection was also obtained 
by A. Edge and P. Salome at the same velocity as the C0(2-l) line, 
suggesting that the CO(l-O) data is beam diluted. This means that 
the CO is probably concentrated in a region smaller than ll"(17 
kpc). The CO line ratio (if correct) is close to or larger than 4, sug- 
gesting this as well. [O'Dea et al.| jl994a| ) detected an unresolved 
narrow H I absorption component in the nucleus with a FWHM of 
~ 220 ± 10 km sec"', as well as a broader component spatially re- 
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Figure 3. //5r imaging of FUV, optical, NIR, and line emission associated with the ~ 10 kpc-scale nebula at the center of the A2597 BCG. The FOV of each 
figure is approximately 10" X 10" 15 kpcx 15 kpc). Panel (a) is FUV continuum emission, attributed to ongoing star formation, from the ACS SBC F150LP 
observation of|Oonk et aLljlOll) . Lya emission is not included in the bandpass. Panel (b) is S-band optical continuum, Ha+[N II], and [S II] emission from 
[Holtzman et al.|jl996^ . Panel (c), also from |Holtzman et al.|jl996") , contains blue optical continuum and a small contribution from [O II]A3727 A emission, 
which dominates the bandpass in the F410M image shown in panel (d), from |Koekemoer et al.jj 199*9) . Panel (e) is primarily emission from the 1.956 fim (rest 
frame) 1-0 S(3) H2 line, originally published by |Donahue et al.lpOOO) . Panel (/) is //-band NIR stellar continuum emission, also from |Donahue et ar]j2000} . 



solved on the scale of the 8.4 GHz radio lobes with a FWHIVI of 
~ 412 ± 40 km sec"' , consistent with the Herschel results. 

[Donahue et aT] pOOO ^ detected H2 1-0 S(3) 2.2 emis- 
sion stemming from vibrationally excited molecular hydrogen on 
roughly the same spatial scale, with some resolved morphology 
tracing the optical emission line filaments. The emission is far too 
bright to be directly accounted for by a cooling flow, and UV irradi- 
ation by the young stellar component was the only model capable 
of accounting for both the HQ/H2 line ratios, as well as the FUV 
continuum strength. That study ruled out AGN photoionization, 
fast and slow shocks, and X-ray heating as the dominant ionization 
source. Excitation by non-thermal electrons may play an important 
role ( |Ferland et al.|2009> . [JaffeetaL ( 2005 1 and Oo nk et al.lpOTO) 
have detected molecular gas out to a maximal radius of ~ 20 kpc 
northward of the nucleus. Thus far, this is the largest known radial 
extent at which a tracer of the warm/cold ISM has been detected in 
A2597 (the longest CC BCG optical filaments have been detected 
out to ~ 50 kpc, for Perseus/NGC 1275, e.g., |Conselice et al.|200T| 
|Fabian|200"3| [Hatch et al.|2006^ . This does not imply that there is 
no warm/cold gas beyond this boundary, as all of these observa- 
tions are limited by sensitivity. The average velocity profiles for 
the ionized and molecular gas, studied by |Oonk et al.| ( (20T0} using 
ESQ Very Large Telescope (VLT) SINFONI data (see Table 1), are 
comparable to one another and are centered at roughly the systemic 
velocity. The data show smaller scale velocity (and velocity disper- 



sion) excesses that appear to be associated with the propagation of 
the radio source ( [Tremblay et al.|2012a^ . 

Gas cooling from an ambient X-ray atmosphere should be 
associated with low net angular momentum. Its velocity profile 
should therefore smoothly increase from the systemic velocity at 
larger radii (e.g., ^ 50 kpc) to a few hundred km sec"' in the center 
of the BCG. As discussed above, this is consistent with what has 
so far been observed. On the other hand, cold gas acquired through 
a merger may be associated with higher net angular momentum 
and exhibit bulk rotation at larger radii. The relative consistency of 
the [O I], C0(2-l), H I, and H2 velocity profiles seems consistent 
with what would be expected if residual ICIVI cooling deposited 
these gaseous components into the BCG with low excess angular 
momentum or velocity structure. While we cannot speak to the ve- 
locity structure of the X-ray gas on these scales, the fact that the 
X-ray isophotes are consistently elongated along the major axis of 
the BCG (top right panel of Fig.[TJ supports the notion that the two 
components are roughly similar in terms of angular momentum. 

Alternatively, a merger origin for the gas cannot be ruled out. 
Although high velocity gas is not detected at large radii, and al- 
though the cold gas is concentrated in the very inner regions of 
the nucleus without complex or high amplitude velocity structures, 
one can imagine a merger producing similar results if enough time 
has passed since the merger for the gas to reach the nucleus and 
couple with the local dynamics. A possible merger origin for the 
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Figure 4. Spectral energy distribution of A2597, from radio (left side) through optical (right side) frequencies. References for the data points used in this plot 
can be found in Table[T] Following the method from |Mittal et al. | j2011^ , a two-temperature blackbody has been fit to the FIR portion of the SED, sampled by 
the new Herschel PACS and SPIRE observations. See Section |4.3| fora discussion of this figure. The 500 /xm SPIRE point is an upper limit. Because the plot 
is log-log, the error bars on most points are smaller than the points themselves. 



warm/cold gas phases in A2597 is discussed in detail by |Koeke^ 
|moer et al.| ( [l999| l. 

A purely merger-based origin for this gas would be difficult 
to reconcile with the FUSE and XMM-Newton results consistent 
with cooling gas < 10^ K, unless conductive interfaces produce the 
lines attributed to cooling (e.g., [Sparks et al. 2012). The observed 
star formation X-ray cooling time threshold (Rafferty et al. 12008) 
|Cavagnolo et al.|2008^ , which we will discuss in Section 5, would 
also be difficult to understand in a "mergers only" scenario. In the 
sections below we will argue that ICM contributions likely play a 
very significant role in A2597, regardless of the unknown contri- 
bution by cold gas acquisition scenarios like mergers. We reiterate 
that a recent gas-rich major merger is almost certainly ruled out for 
A2597. 



4.3 New Herschel constraints on the cold gas and dust 
components 

The 1.8 ± 0.3 x lO' Mq cold molecular gas reservoir in A2597 is 
cospatial with the ~ 10 kpc emission line nebula, which has been 
studied extensively in the literature (e.g., |Heckman et al.|I989[[Mc-| 



Namara & 0'Connell"1993'; 'Voit & Donahue" 1997" "McNamara 
et al..,1999,.Koekemoer et al. J999, .Oegerle et al., 2001, .O'Dea 



et al.|2004[paffe et al.|2005[|Oonk et al.|2011| l. High spatial reso- 
lution archival HST imaging, shown in Fig. [3] reveals its complex 
filamentary morphology in the FUV and optical (panels a through 
d in Fig.|3j. Deep optical spectroscopy by | Voit & Donahu^jI997| > 
shows the gas temperature, abundance, and electron density of the 
nebula to be 9,000-12,000 K, ~ O.SZq, and n, ~ 200 cmr\ respec- 
tively. [Koekemoer et al.| ( |1999^ estimated the pressure associated 
with the 10" K gas on these scales to be 



Pnebula ~ 5.5 X 10 



200cm-V VlO* K 



dyn cm 



(9) 



Our estimated central X-ray pressure on the s ame ~ 10 kpc scale 
is 7 X 10"'" dyn cm"^ (Tremblay et al. |2012a| , very close to what 
is inferred for the lO" K gas, implying that the two phases are in 
rough pressure equilibrium. 

The nebula is associated with a substantial dust component 
that can be quantitatively studied using the Herschel observations. 
In Fig.|4]we show the radio through optical SED of A2597, with 
the FIR component filled in by the new Herschel PACS and SPIRE 
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data (green circles). Following the same method described in |Mit-| 
|tal et al.| ( [201 \) , a two-temperature black body is fit to the Herschel 
(PACS and SPIRE) and Spitzer 24 jim data points to model the con- 
tribution by dust to the FIR SED. Our fit assumes a dust absorption 
coefficient of = 5. 6 x (i//3000 GHz) ^ m^ kg"' with a dust emis- 
sivity index of /? = 2 ( |Dunne et al.|200o[|Mittal et al.|201 1) . We note 
that the 850 /im SCUBA data point could also be associated with an 
additional cold dust component, though we choose not to include 
it in our modeling given the uncertain radio source contribution to 
this flux. 

We find the best-fit temperature and mass of the warm 
dust component to be Twarradust = 47 ± 1.4 K and M„aimdust = 
(1.7±0.6) X 10^ Mq, respectively. The temperature and mass of 
the cold dust component is found to be Tcoid dust = 20 ± 1.7 K and 
McoW dust = (1.3 ±0.5) X 10^ Mq, respectively. This two component 
blackbody fit could be unrealistic if (1) there is a large population of 
small dust grains driving a broad emission distribution that might 
mimic a modified blackbody, and/or (2) if the warm dust grains 
are heated stochastically (which is only realistic for very small 
grains, e.g., |Draine & **2001). Furthermore, a two-temperature 
fit assumes that there is no dust at temperatures between ~ 20 K 
and ~ 48 K, which is almost certainly not true. The SPIRE and 
SCUBA points suggest the presence of very cold dust, which would 
increase the mass limit. Finally, the best-fit dust mass is extremely 
sensitive to the assumed dust temperature, which in turn is degen- 
erate with the dust emissivity index j3. We have assumed that /3 = 2, 
but the value can vary between ~ 1.5 and 2 (e.g., |Tabatabaei et al.| 
|201 l i While we find that varying /3 between these values does 
not significantly change our results, we stress that the dust masses 
and temperatures listed here are extremely assumption-heavy. Al- 
tering these assumptions even slightly can significantly alter the 
result. As a point of illustration, the single-temperature fit to the 
A2597 Herschel points by |Rawle et"S^pOT2l ) finds T" = 31 K and 
Mdust = 1.6 X 10** Mq, an order of magnitude higher than our re- 
sult. Ultimately, we lack the data to speak affirmatively about the 
validity of these many assumptions. 

With these disclaimers noted, our total Herschel-deriwed dust 
mass of Mcoid dust = 1.3 x 10^ Mq compared to the inferred molec- 
ular gas mass of 1.8 ± 0.3 x lO' Mq implies a gas-to-dust ratio of 
140. We again stress that, while Herschel can place the tightest 
available constraints on the A2597 dust mass, this gas-to-dust ratio 
remains an uncertain lower-limit given (1) the dust mass, tempera- 
ture, and fit parameter degeneracies discussed above, and (2) CO- 
inferred molecular gas masses are inherently uncertain because of 
the CO-H2 conversion "X factor" (e.g., [Liszt et al.|2010[ and refer- 
ences therein). This is especially true considering the apparent low 
metallicity of the gas, which could mean that our assumed CO-H2 
conversion factor is an order of magnitude too low (e.g., |Bolatto| 
|et al.|[20TT) . We could therefore be significantly underestimating 
the gas mass. 

Nevertheless, the inferred lower-limit ratio is comparable to 
Galactic ratios, and somewhat low relative to past estimates for 
other CC BCGs (e.g., |Edge,2001 j . These older ratios suffer from 
an inability (at the time) to place convincing limits on CC BCG 
dust masses, so a more comprehensive study with the new Her- 
schel sample is needed to better understand the "typical" gas-to- 
dust ratios of CC BCGs (assuming such a typical ratio exists — it 
certainly may not). Unpublished and science demonstration phase 
results for the Herschel CC BCG sample so far indicate more MW- 
type gas-to-dust ratios on average (e.g., |Edge et al.|2010b|a|[Rawle| 
|et al.|2012[ Oonk et al. 2012 in prep). 



The Herschel [C Il]-to-FIR luminosity ratio, which is a tracer 
of the relative cold gas and dust cooling rates, is ~ 0.014 for A2597. 
This value is high relative to similarly FIR luminous star forming 
galaxies (whose [C II]/Lfir ratios are generally an order of mag- 
nitude lower). The ratio seems to be higher for lower metallicity 
galaxies (e.g., |Maiolino et al.|2009| and references therein), so the 
relatively high ratio estimated for A2597 could be partially due to 
lower metallicity gas, consistent with a cooling flow scenario. 

The MW-type gas-to-dust ratio means that the star forming gas 
in A2597 is quite dusty. This is inconsistent with a scenario wherein 
already dusty gas originates in the hot atmosphere, where grain 
sputtering timescales are short (e.g., |Draine & Salpeter|[T979[ re- 
view by Draine 2003 ). In this case, one might expect far higher gas- 
to-dust ratios (i.e., dust-poor gas) stemming from cooling flow ori- 
gin scenarios. Furthermore, even if the cooling flow gas were dusty 
despite short grain destruction times, then extended FIR emission 
might be expected on the same scales as the X-ray gas, as dust 
emission should be a major channel of energy loss for the cooling 
ICM. This extended FIR emission cannot be detected in A2597 by 
Herschel due to spatial resolution limitations ( [Edge et al.|2010a[b| ). 

The above result suggests that while the cooling ICM may 
contribute the bulk of the cold gas in A2597, it cannot contribute 
the substantial dust component associated with that cold gas. In- 
stead it is more likely that the bulk of the dust component is con- 
tributed by stellar mass loss (e.g., via dust-rich AGB winds and su- 
pernovae). Modem cooling flow models suggest that cooling from 
the ICM concentrates in thermally unstable clouds and filaments 
(e.g., Fabian et al. 2011) . In these thermal instabilities, the cool- 
ing/cold gas may shield the dust from sputtering. The recent detec- 
tion by [Donahue etlT] pOTT) of PAH emission in A2597 almost 
certainly necessitates such a shielding mechanism, given the ex- 
treme fragility of PAH molecules in hot ambient environments. If 
the dust is not efficiently shielded within the filaments, then dust 
production rates must be higher than what is implied by most mod- 
els for dust production. 

[Voit & Donahue[|201 1^ argue that mass loss from the old stel- 
lar population of the BCG is indeed a significant source of dusty 
gas in many cool core clusters. We note that this scenario is not 
at odds with the residual cooling flow model for A2597. As fila- 
ments condense from the cooling hot atmosphere, they do so within 
the stellar body of the galaxy and should therefore contain existing 
stars within their volume. [Voit & Donahue [ ( (20 1 Ij ) suggest that the 
ejected dust-rich envelopes stemming from the old stellar compo- 
nent are not assimilated into the hot phase, but rather remain cold 
and confined to the cold ambient clouds and filaments, where the 
dust may be shielded from interaction with the hot phase. If the old 
stellar component plays an important role in enriching the filaments 
with dust, then it is likely that it at least plays a non-negligible role 
in contributing to the observed cold gas mass as well. 



5 NEW RESULTS ON STAR FORMATION 

5.1 Herschel constraints on the star formation rate 

The 8-1000 /im luminosity can be used with the Kennicutt relation 
to estimate the FIR-inferred star formation rate l |Kennicutt| I998| l; 

;i|^s«x(,5^), uo, 

Using Lfir = (6.50 ± 1.35) X 10*' ergs s"' (1.69±0.35x lO'" Lq) 
derived from the Herschel data, we find SFRfir ^ 2.93 ±0.69 Mq 
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Figure 5. (left) X-ray cooling time profile vs. cluster-centric radius. The green line marks the ?Ti/fff <^ 10 thermal instability threshold discussed by Sharma 
|et al .ll 201 1 1, where ?ti is the thermal instability timescale and % is the local gravitational freefall timescale. [Sharma et al.H2011^ suggests that, beyond this 
threshold, a cooling flow can form thermally unstable, star forming clouds and filaments. The blue shaded region marks the largest radial extent of observed 
ionized and molecular gas. The red shaded region, which overlaps it, marks the largest radial extent at which X-ray cavities are observed, {right) Ratio of the 
thermal instability timescale and freefall timescales vs. cluster-centric radius. The same fTi/fff <^ 10 star formation threshold is marked with a green line. See 
SectionB^for a discussion of these profiles. 



yr '. This is consistent with the many past SFRs estimated for 
A2597, which range from 2 - 12 yr"' depending on the method 



used jMcNamara & 0'Connell|1993| [O'Dea et al.|2004i |Donahue| 
let al.|2007[|Oonk et al.|20TT]|Rawle et al.|2012[ l. Estimates of star 
formation rates (particularly in the FTJV) are extremely sensitive 
to internal extinction by dust. The Balmer sequence in A2597 sug- 
gests that the internal extinction is significant (about one magnitude 
in V-band, i.e., Av ~ l, |Voit & Donahue|1997||Oonk et al.|20TT] l, 
consistent with the substantial dust mass inferred from the Herschel 
data. 



5.2 Star formation entropy threshold 

Using numerical simulations and ignoring thermal conduction, 
ISharma et al.| ('20Hi showed that local thermal instabilities in a 
cooling flow will only produce a multiphase, star forming ISM 
when the ratio of the thermal instability timescale fxi to the lo- 
cal gravitational free-fall timescale fff is fxi/^ff ^10 (see 



Sharma 

|et al.|(20TT| > for definitions of these quantities). Expressed m terms 
of the gas entropy S, this thermal instability threshold is 5 ~ 20 
keV cm^. Importantly, this theoretical result is very close to the ob- 
served 5 ~ 30 keV cm^ star formation onset threshold discussed by 
[Rafferty et a"Ll ( |2008) ; [Cavagnolo et afjpOOS) . 

The fxi/fff ^10 threshold is marked in the green line on the 
X-ray cooling time profile in Fig. [SJj as well in Fig. [5^, in which 
we plot the fri/fff ratio vs. cluster-centric radius. In both panels we 
mark maximal radial extent of (a) Lya, FUV continuum associated 
with star formation, as well as warm and cold molecular gas in- 
ferred from Very Large Telescope (VLT), HST, CO and Herschel 
observations and (b) the X-ray cavity network with blue and red 
shaded regions, respectively. Note that the X-ray cavity network 
(red region) overlaps with the blue region, meaning it also extends 



inward to the cluster centre (as can be seen in the bottom two pan- 
els of Fig. [TJ. We calculate the fxi/te ratio from the X-ray data 
by fitting third-order polynomials in log space to the temperature 
and pressure profiles, then analytically differentiating the logarith- 
mic pressure profile to obtain the free-fall time. We account for the 
presence of the BCG by setting a minimum value for the gravita- 
tional acceleration g to be that of a singular isothermal sphere with 
a velocity dispersion of 250 km s"' . This BCG contribution is only 
important at radii 10 kpc. 



As discussed in [Sharma et al.|(2011| l, the thermal instability 
timescale can be approximated as a multiple of the cooling time 
that is dependent upon how much AGN heating one includes in 
the model. While the curve in Fig.|5j? is for pure free-free cooling, 
including a heating component could scale it upwards by a factor 
of ~ 2 (for moderate heating) to ~ 3 (for strong heating). Regard- 
less of whether or not heating is included, the [Sharma et al.[ l |201 l| l 
predictions are roughly consistent with observations in the case of 
A2597. The same is true of the 30 keV cm"^ star formation entropy 
threshold discussed by [Rafferty et al.[ l |2008^ and jCavagnolo et aT] 
POOS^ — we don't mark this threshold, but it would lie in nearly 
the same locations as the green lines in Figs.|5j? and b. Therefore in 
either case, the threshold lies just outside the region containing star 
formation, ionized and molecular gas, and X-ray cavities. 

This apparent threshold would be difficult to understand if 
the cold gas fueling star formation was supplied by a merger. One 
might imagine such a scenario if there were a steep radial gradi- 
ent in thermal conduction efficiency, such that clouds introduced 
by a merger no longer evaporate at the radius where conduction no 
longer outpaces radiative cooling. This is purely speculative, how- 
ever, and the results presented here strongly suggest that a residual 
cooling flow from the ambient hot atmosphere is the primary con- 
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Figure 6. (top left) In red/orange we show the //ST/STIS FUV observation of extended Lya emission associated with the A2597 BCG, from |0'Dea et al.] 
J2004|. 1.3 GHz, 8.4 GHz, and 330 MHz radio contours are overlaid in blue, white, and black (dashed), respectively, while adaptively smoothed 0.5-7 keV 
X-ray contours are overlaid in green. We have removed the innermost 1.3 GHz, 330 MHz, and X-ray contours to aid viewing, (top right) HST/ ACS SBC 
F150LP ~ 8 ksec exposure of FUV continuum emission associated with the A2597 emission line nebula. In red contours, we outline the boundary within 
which low surface brightness FUV continuum emission has been detected at ~ 2a above the background, (bottom left) The unsmoothed 150 ksec 0.5-7 keV 
X-ray data, shown with an aggressive color scale stretch used to highlight the high surface brightness knots observed near the X-ray centroid. Lya contours are 
overlaid in red. The northern high surface brightness X-ray knot bends upwards to follow the overall distribution of the high surface brightness Lya and FUV 
continuum emission, (bottom right) The same X-ray image, with red boxes mai'king the sectors from which X-ray spectral data were extracted and modeled. 
A WABS X MKCFLOW fit to the data finds a combined mass deposition rate of ~ 7 it 0.5 Mq yr"' associated with the X-ray filaments. This is comparable to the 
local star formation rate estimated from FUV emission in the filaments. 



tributor of the warm and cold gas in A2597. We provide further 
supporting evidence for this suggestion in the following section. 



5.3 Evidence for filamentary cooling cliannels 

In Fig. |6|3 we show the extended filamentary Lya emission 
(HST/SJIS, from |0'Dea et al.|2004^ , whose brightness can be ac- 



counted for by photoionization from the underlying FUV contin- 
uum emission stemming from the young stellar component, shown 
in Fig.|6}7 (HST/ACS, from Oonk et al. 2010). In the left panel, we 
show the adaptively smoothed 0.5-7 keV X-ray contours in green 
(with the central contours removed to aid viewing), and the 330 
MHz, 1.3 GHz, and 8.4 GHz radio contours in black, blue, and 
white, respectively. In the right panel we use red contours to mark 
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Figure 7. (left) Unsharp mask of the 8 ksec HST/ ACS FUV continuum image shown in Fig.|6] Note the strong spatial correspondence of the Northern arc of 
FUV emission with the leading edge of the 8.4 GHz radio lobe, which is outlined in red contours. There is also evidence that the southern radio lobe has swept 
out the star forming gas, or has been allowed to expand after emerging from the dense gaseous medium, (right) Single Stellar Population (SSP) age map from 
|Oonk et ah 1 201 1 1, made by comparing observed FUV/J/-band ratios (from HST and the VLT, respectively) to those predicted from Brazual & Chariot 1 20(3} 
models. The youngest portion of the young stellar component is found nearer to the central 8.4 GHz radio source, which we overlay in white contours. The 
young stars also extend Northwards along the region threaded by hya and FUV continuum filaments (as can be seen by comparing this figure with Fig.|6}. 



the outermost boundary at which low surface brightness FUV con- 
tinuum emission is observed at a level of 2a above the background. 

The faint underlying FUV continuum follows the extended 
bright Lya clumps and filaments northwards and southwards, out 
to a cluster-centric radius of nearly 20 kpc. In the X-ray tempera- 
ture map presented in Tremblay et al. ( 2012a l, there is an elongated 
region of cold X-ray gas extending along the same N-S projected 
axis. This filamentary axis is also significantly offset in position 
angle from the cavity/radio "AGN heating axis". The Northeastern 
Lya filaments lie along a rim bordering the cold X-ray filament 
(feature 5 on Fig. [TJ, which could be due to the star forming gas 
having been swept outward by the 1.3 GHz radio source (or by 
the soft X-ray filament that the radio source may have dredged up- 
wards, e.g., [Tremblay et al.|2012a| >. 

In the bottom two panels of Fig.|6]we show the 0.5-7 keV un- 
smoothed X-ray data, with an aggressive color scale stretch applied 
to enhance the contrast. Lya contours are overlaid on the bottom 
left panel. Note how the Northern bright X-ray knot bends 5 kpc 
upwards to follow the high surface brightness Lya and FUV con- 
tinuum emission. The red boxes in the bottom right panel of Fig. [6] 
right panel mark the sectors from which we extract and model X- 
ray spectral data. These sectors approximately cover the region over 
which the extended Lya and FUV continuum filaments are cospa- 
tial with the soft X-ray gas. As stated previously, this is also the 
axis along which elongated cooler X-ray gas is observed in the X- 
ray temperature map. Following the method described by |Tremblay| 
|et al.| ( |2012a) l, a cooling flow model (WABS XMKCFLOW) fit to the 
spectra extracted from these sectors finds a combined mass deposi- 
tion rate of ~ 7 ± 0.5 Mq yr"' associated with the region of X-ray 
gas cospatial with the ionized filaments filaments. This is roughly 
comparable to the local star formation rate estimated from FUV 



emission in the filaments. A multi-component fit to the X-ray spec- 
tral data (e.g., S anders et al.||20(j4} suggests that multiphase gas 
could be present, though the quality of the fits prevents any con- 
clusions from being made. That these regions could account for a 
significant fraction of the residual cooling luminosity is consistent 
with recent similar results on Perseus and Abell 2146 from lFabianI 
|etal.| l |201l) and Russell et al. (2012), respectively. 

This result serves as further supporting evidence to our overall 
suggestion that (a) ICM cooling contributes a significant fraction of 
the cold gas budget fueling the star formation and (b) this cooling 
happens in a spatially discrete, structured manner concentrated in 
filaments. In the case of A2597, these filaments may preferentially 
be found in a cooling channel that is nearly perpendicular to the 
projected AGN heating axis permeated by X-ray cavities and radio 
emission. Future Atacama Large Millimeter/submillimeter Array 
(ALMA) observations capable of spatially resolving molecular gas 
emission on the scale of the filaments will be necessary to either 
refute or reinforce the above claims. 



5.4 Evidence for both triggered and persistent star formation 
amid AGN feedback 

The frequently studied morphological correspondence of blue ex- 
cess filaments with the 8.4 GHz radio source in A2597 presents 
an important test case for models of jet-induced star formation. 
The morphological evidence can be seen in the leftmost panel of 
Fig.|7] in which we show an unsharp mask of the FUV continuum 
emission. The leading edge of the northern 8.4 GHz radio lobe, 
shown in red contours, spatially correlates with the northernmost 
FUV arc, which might be expected if star formation is triggered by 
shock-induced cloud collapse as the propagating plasma entrains 
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and displaces cold gas phases (see e.g., the shock/jet-induced star 
formation m odels by|Elmegreen & Elmeg reen|1978[|Voit|1988|[De| 
|Young|1989||McNarnara & OX :onnell 1993). 

The optical IFU spectroscopy of |Oonk et al.|pOTO) may lend 
evidence in support of this, as it shows high velocity dispersion 
molecular gas along the "sweep-up" trajectory leading to the north- 
em FUV arc. Further evidence for gas displacement by (or per- 
haps dynamical confinement of) the radio source can be seen in the 
southern lobe, which anti-correlates with the local FUV distribu- 
tion. It is also worth noting that the "U" shaped loop of filamentary 
star formation directly southwest of the southern radio lobe is ex- 
tended along the general axis of the 330 MHz radio source and the 
large western X-ray cavity. The compelling morphological corre- 
spondence has motivated several past investigations of jet-triggered 
star formation in A2597 i jDe Young|1995[|Koekemoer et al.|1999| 
|0'Dea et al.|2004) . Each of these works found the model energeti- 
cally feasible on timescales that make sense in relation to the radio 
source ( |De Young|19 95 1. 

It is worth pursuing these timescale estimates further, by com- 
paring the estimated radio source and X-ray cavity age limits 
with the young stellar component dating presented by |Oonk et al.| 
( |20II) , using data from HST&nA the VLT FORS imager. That work 
compared measured FUV/C -band ratios with those predicted from 
|Bruzual & Chariot] j2003| l simple stellar population models. We 
show the YSP ages inferred by jOonk et al.| ( [20TT| l in the rightmost 
panel of Fig. |7] Note that these estimated ages are better consid- 
ered in the relative rather than absolute sense. The ages are un- 
avoidably uncertain given the age-metallicity-extinction degener- 
acy associated with measuring ages of a young stellar component. 
The FUV/(/ -band color is strongly dependent upon the nature and 
patchiness of the intrinsic reddening, which cannot be quantified on 
these scales. Furthermore, while C IV A1549 A resonant line emis- 
sion has so far not been convincingly detected in A2597 ( |0'Dea| 
|et al.|2004{ |Oonk et al.|20Il^ , a significant contribution from this 
line would result in an underestimation of stellar ages from the 
FUV/C-band color. We note that C IV has recently been detected 
in the //5r Cosmic Origins Spectrograph (COS) FUV spectroscopy 
of M87 ( [Sparks et al.|20T2| >. 

With these caveats in mind, the age map shows younger stars 
nearer to the radio source. Looking on smaller scales, there is possi- 
ble evidence of even younger (5 Myr) stars found (a) at the northern 
edge of the northern radio lobe (where jet-triggered star formation 
could have recently occurred), and (b) along the projected axis of 
the host galaxy stellar isophotes. When we compare X-ray cavity 
ages from this paper with these YSP ages, we find that the range 
of young stellar ages entirely encompasses the inferred age range 
of the X-ray cavity network. |Oonk et aTlpOI l| l found ages for the 
YSP ranging from 5 - 700 Myr old, while our estimates for the X- 
ray cavity ages range from 10-200 Myr. The stellar ages also en- 
compass the age range estimated for the radio sources (see Section 
|3.4| l. The stellar ages inf erred by|Oonk et al.|pOTT] l are consistent 
with similar estimates by |Koekemoer et aL ( 1999^ . 

Two independent studies infer a young star age range with 
limits that are older and younger than the age limits for the X-ray 
cavity network. While acknowledging the many caveats associated 
with these estimates, this could imply that low levels of star for- 
mation (2- 12 Mq yr"') have managed to persist even amid the 
AGN feedback-driven excavation of the X-ray cavities. The stellar 
age spatial distribution may further suggest that star formation was 
more spatially extended in the past. This interpretation is ambigu- 
ously dependent upon the likely fact that most star formation occurs 
in the filaments, which are more concentrated in the center than in 



the periphery. Nevertheless, we cannot rule out the possibility that 
even if the star formation has been persistent during the current 
AGN feedback episode, it may be qualitatively changing with time 
(i.e., declining and becoming more spatially concentrated). 



6 CONCLUDING DISCUSSION 

The results presented in this paper, considered in the context of 
three decades of previous work on A2597, motivate the following 
general conclusion: while radio-mode AGN feedback has injected 
enough energy into the hot ICM to inhibit the classical cooling flow, 
the source harbors a residual cooling flow at 4% to 8% of the ex- 
pected classical rates, giving rise to star formation amid the ~ lO' 
Mq cold molecular gas reservoir in the nucleus. While we cannot 
rule out cold gas contiibutions from mergers or tidal stripping, we 
suggest that the cooling ICM is likely the dominant supply channel 
for cold gas to the nucleus in A2597. 

A simple test of this model can be made by considering the 
predicted residual cooling flow mass deposition rate with the ob- 
served mass of the cold molecular gas reservoir, along with the 
star formation and black hole accretion rates. Together, these can 
be considered the ultimate mass sinks of the residual cooling flow, 
so their mass and energy budgets should be consistent with one 
another. In Table 3 we compile the various mass, mass flux, en- 
ergy, and timescale limits assembled throughout the course of this 
analysis. We separate out the important quantities related to inhibi- 
tion of the classical cooling flow (i.e., the classical mass deposition 
rate and the lower limit kinetic energy input estimate based on X- 
ray cavity analysis). We then list the multiphase constraints on the 
residual cooling flow model. The predicted residual mass deposi- 
tion rates are roughly 20-40 Mq yr"' within 30 kpc. A steady 
residual mass deposition rate of this magnitude would accumulate 
the observed ~ lO' Mq of gas in the central 30 kpc within <_ 10* 
yr. Taking the ratio of the observed molecular gas mass to the ob- 
served star formation rate yields a gas depletion timescale that is 
also on the order of 10* yr. It is therefore possible for there to be an 
approximately steady-state distribution of multiphase gas masses 
over the AGN lifetime, with a possible slow accumulation of cold 
gas at a rate comparable to the difference between the local mass 
deposition and star formation rates. 

While the proposed model appears successful in this regard, 
the above argument is of course very qualitative. [Gaspari et al.| 



l |2011i presented a follow up study to the Sharma et al. 1 201 1 



suits, predicting (among other things) the mass budgets expected 
from various cooling flow recipes, over time, in the innermost 20 
kpc of a BCG. Their results for an uninhibited cooling flow expect- 
edly overpredict the mass of the warm and cold phases by several 
orders of magnitude. When they include AGN heating at varying 
efficiencies, their results fall more in line — albeit on the high end 
— with what is observed in A2597. Their feedback-inhibited resid- 
ual cooling flow models predict a roughly steady accretion rate of 
~ 10 Mq yr"' of r < 5 X 10^ K gas within the central 20 kpc. Af- 
ter 1 Gyr, this accretion rate would accumulate a ~ 5 x lO'" Mq 
cold gas reservoir, which is an order of magnitude more massive 
than what has been measured in A2597. However, as |Gaspari et al.| 
(j20TTjl notes, their simulations do not include star formation, which 
could reduce the accumulation rate of cold gas by an order of mag- 
nitude. On roughly the same scale, the highest measured SFR for 
A2597 is ~ 10 M© yr"' , which is roughly the cold gas accretion rate 
predicted on the same scale by |Gaspari et al.| l [2011^ . It is also simi- 
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Table 3. A summary of the various mass, mass flux, energy, and timescale limits compiled as part of this paper. All values are measured or 
estimated at or within a central radius of 30 kpc. (1) Temperature range of the ISM phase (if applicable); (2) qualitative description of that 
phase or feature (3) mass associated with the phase or feature; (4) power associated with the phase or feature; (5) estimated mass conversion 
flux or deposition rate; (6) A rough timescale associated with the listed phase or feature, (cool is a cooling time associated with the listed 
cooling flow, fheat is a crude estimate of the possible timescale over which X-ray cavity enthalpy could be dissipated as heat in the ISM, 
limited by the cavity lifetime, (deplete is the cold gas depletion timescale, set by the ratio of the observed cold molecular gas mass and the 
star formation rate, faccum is the time it would take the predicted residual cooling flow mass deposition rate to accumulate the observed cold 
molecular gas mass. Note that the coiresponding timescales (t^aol vs. fhcat and 'accum vs. fdcpiete) roughly balance one another. 
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A/phase(K^30kpc) L(R^30kpc) 


M (r ^ 30 kpc) 
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(Mq) (ergs sec"') 
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(2) 


(3) (4) 


(5) 


(6) 




Limits on a Classical Cooling Flow 


10' -lO** 


Classical Cooling flow 


(l-4)xl0'"' 
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30 




AGN Feedback Input 


> 1.89x10"** 






-40 


Limits on an Residual Cooling Flow 


10' -lO* 


Residual Cooling flow 


(0.4-3.2) X 10« 


Meool ~ 20-40 


^accum 


- 10 


10« 


FUV Cooling flow 


^4x10"° 
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lO'^-lO^ 


Warm Ionized gas 


> 9.7±0.3x lO^ > 3.5x10"! 


MsFR~2-12 






10- 10^ 


Warm/cold molecular gas ^ 1 .8 ± 0.3 X lO' ^5x10"' 




^deplete ^ 


- 10 



lar to the X-ray derived mass deposition rate for the region cospatial 
with the filaments of ionized gas, as discussed in Section[53] 



7 SUMMARY 

We have presented a multiwavelength study of the central brightest 
cluster galaxy in the cool core cluster Abell 2597. The main results 
of this paper can be summarized as follows. 

• New Chandra observations reveal the X-ray cavity network 
to be more extensive than previously known, and associated with 
enough enthalpy to locally inhibit the classical cooling flow. 

• A comparison of estimated cavity and radio source ages sug- 
gests that the AGN duty cycle is near to 100%, requiring a near- 
steady sub-Eddington flow of gas to the nucleus. 

• The //ericAeZ-derived warm dust temperature and mass is esti- 
mated to be r„arm dust = 47 ± 1 .4 K and M„arm dtsl = (1 .7 ± 0.6) x 10^ 

Mq, respectively. The temperature and mass of the cold dust 
component is found to be Tcoid du,st = 20 ± 1.7 K and Mcoid dust = 
(1.3 ±0.5) x 10^ Mq, respectively. We discuss important uncer- 
tainties and assumptions associated with these values in Section 
4.3. 

• We present an updated CO-inferred cold molecular gas mass 
of 1.8 ±0.3 X lO' Mq, using previously unpublished CO(2-l) 
IRAM 30 m observations. 

• The newly measured gas and dust masses yield a relatively 
low. Galactic-type gas-to-dust ratio of ^ 140, so the gas is dusty. 
We argue that mass loss from evolved stars is the most likely source 
of this substantial dust component. A shielding mechanism is likely 
required to protect the grains from interaction with the ambient hot 
ISM. 

• The warm and cold gas phases are approximately cospatial 
with the coolest regions of X-ray gas. 

• Dynamical constraints on the central cold gas component do 
not permit large scale rotation or major asymmetric velocity struc- 
tures, consistent with a scenario wherein the cold gas is accreted 
with low net angular momentum, inconsistent with what is ex- 
pected from a merger. 



• The theoretical thermal instability threshold lies just outside 
the observed <; 30 kpc X-ray cavity network. Molecular and ion- 
ized, star forming gas lies interior to both. This result is consis- 
tent with theoretical predictions of the entropy threshold at which 
ICM cooling begins to form thermally unstable cold clouds and fil- 
aments. 

• X-ray derived mass deposition rates along regions cospatial 
with star forming filaments are consistent with the locally estimated 
star formation rates, strongly suggesting a causal connection. 

• The young stellar component occupies an age range that is ap- 
parently wider than that for the X-ray cavities. This could suggest 
that star formation has persisted amid the feedback-driven excava- 
tion of the X-ray cavity network. Localized sites of star formation 
may also have been triggered by the propagating radio source. 

We conclude that a residual cooling flow with a strength of 4% 
to 8% of the expected classical mass deposition rates is the domi- 
nant contributor of the cold gas reservoir fueling star formation and 
AGN activity in the AbeU 2597 BCG. 
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